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Abstract 

Within a simple quantization scheme, observables for a large class of finite 
dimensional time reparametrization invariant systems may be constructed 
by integration over the manifold of time labels. This procedure is shown 
to produce a complete set of densely defined operators on a physical Hilbert 
space for which an inner product is identified and to provide reasonable results 
for simple test cases. Furthermore, many of these observables have a clear 
interpretation in the classical limit and we use this to demonstrate that, for a 
class of minisuperspace models including LRS Bianchi IX and the Kantowski- 
Sachs model this quantization agrees with classical physics in predicting that 
such spacetimes recollapse. 
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I. INTRODUCTION 



A generalization of Kuchaf's "Hilbert space problem" [ l|] provides one of the greatest 
obstructions to developing a quantum theory of gravity. At issue are the questions "What 
are the physical (quantum) observables?" and "On what Hilbert space (if any) do they 
act?" This includes the familiar question "What is the inner product on the space of phys- 
ical states?" Here, we combine the idea emphasized by DeWitt [[| of constructing "local" 
diffeomorphism invariants by integrating singular densities over a manifold with the quanti- 
zation of constrained systems described by Dirac M and find that we are able to construct 
gauge invariant operators which (at least in the classical limit) have a clear interpretation 
and which are densely defined on a Hilbert space of physical states with a known inner 
product. This approach has the advantage that it is not necessary to identify an "internal 
time function" |IJ, that it may be applied to a large class of minisuperspace models - which 
includes the LRS Bianchi IX and Kantowski-Sachs cosmologies - and that for the models 
just mentioned it predicts that the corresponding spacetimes recollapse, in accord with the 
classical result. Note that it is highly nontrivial for a quantization scheme to satisfy this last 
condition and, in particular, that procedures based on a Klein-Gordon type inner product 
instead describe the universe as expanding "forever" [Q. 

Let us consider for a moment the more general situation of a diffeomorphism invariant 
system on a spacetime manifold M. Some interesting common examples are general rel- 
ativity, topological field theories, parametrized Newtonian systems, and homogeneous cos- 
mological models (such as the Bianchi models) for which M is just the manifold R of time 
labels. Solving the equations of motion for such systems introduces a number of arbitrary 
functions corresponding to the freedom to choose coordinates on M. Physical information is 
contained in gauge invariant quantities constructed from such solutions (that is, quantities 
which do not depend on these arbitrary functions) so that, classically, we face the issue of 
constructing such objects and extracting useful information. We will refer to such gauge 
invariants as observables. 
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As has been pointed out a number of times (see, e.g. 0), we do in fact know how 
to construct a great number of observables for classical diffeomorphism invariant theories. 
Given any density uj on an n-manifold M, J M to is such an observable. Of course, these 
observables are not of the type with which we are familiar from, say, scalar field theory and 
they may be highly non-local. This is the case for simple examples such as f M y/—gR in 
Einstein gravity. Also, if the density u is not chosen carefully, such an integral may converge 
only when evaluated for very special solutions of the dynamics. 

However, as described in f2|, some observables of this type can be easily interpreted. 
If the density u is distributional, the resulting gauge invariant may be effectively local on 
M. A simple example of such an observable is the value of some scalar quantity at a point 
specified by an "intrinsic coordinate system." For gravity coupled to a set of scalar fields <p k , 
k G {0,1,2,3}, the corresponding integral might be 



taneously only at a finite number of spacetime points and for which the metric is smooth 
whenever cf) k = 0. 

We would like to use similar ideas to define and discuss observables of this type for 
quantum theories with reparametrization or diffeomorphism invariance. Two standard crit- 
icisms of observables of this type in the classical setting are 1) that they will in general 
be well-defined (and finite) only on a part of the classical solution space and 2) that, as 
they involve integration over time, they are likely to be of little use unless the theory can 
be solved completely (so that the integrand may be written explicitly). Nevertheless, we 
will see that quantum analogues of these observables exist as densely defined operators on 
a space of physical states for interesting models and that certain properties of the resulting 
operators can be studied even though the models may not be exactly solvable. In par- 
ticular, such observables can be used to show that, in correspondence with the classical 




(1.1) 



where R is the curvature scalar and is a Jacobian. This object defines a well-behaved 
function on any part of the classical space of solutions for which the fields (p k vanish simul- 
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result, our quantization of appropriate minisuperspace models describes spacetimes whose 
homogeneous slices expand to some maximal size and then recollapse. 

The strategy is as follows. First, section [IT] discusses the classical form of the observables 
of interest. This will provide insight into the definitions and arguments used in the quantum 
treatment. We then give the general quantum formalism for the operators and ths physical 



Hilbert space in |TJ and apply it to two test cases in |TV|. After verifying that appropriate 
results are produced for these simple cases, we turn to the cosmological models of interest in 
|V[ In each application, the integrals that define our observables will be shown to converge 
and we find a "complete set" of observables that are densely defined on the physical Hilbert 
space and symmetric with respect to a physical inner product. In addition, [V] shows that 
when applied to "separable semi-bound" cosmological models, this quantization describes 
spacetimes that recollapse. We conclude with a short discussion in section [V]]. 



II. OBSERVABLES THROUGH INTEGRATION 

Section |TJ will address quantization of reparametrization invariant models and will focus 
on the properties of certain observables. In order to provide an intuitive understanding of 
the definition and use of these observables, we now discuss their classical counterparts. Thus, 
this section can be regarded as essentially heuristic background for the more rigorous results 
of[mj |V], and 0. 

We are interested in systems that are time reparametrization invariant; that is, systems 
defined on some d + 1 dimensional spacetime M for which there is a gauge symmetry that 
takes f(t) to any f(T(t)) whenever dT/dt > and / is a scalar on M. We concentrate on 
the case d = and on observables like |1 . 1| . However, because canonical techniques will be 
used, it is better to consider invariants of the form 

[a] q=T = r dt^5(q(t)-r)a(t) (2.1) 

J — oo CLt 

which may be regarded (when the integral converges) as either functions on the space S of 
classical solutions or, through the evolution map that constructs a solution from a piece of 
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initial data specified by a point on the constraint surface, on the part of the corresponding 
phase space in which any constraints are satisfied. When q(t) takes the value r once along a 
solution, this object gives the value of a(t) when q(t) = r up to a sign. Similar observables 
were used in 

The advantage of |2.1| is that it is invariant not just under reparametrizations (t — > T(t) 
with dT/dt > 0), but also under pull backs of the manifold of time labels through maps of 
degree 1 (t — > T(t) for any smooth T(t) of degree 1). This large class of transformations 
typically become gauge symmetries as well when a time reparametrization invariant system 
is written in Hamiltonian form. 

In fact, if the canonical description is in terms of some phase space coordinates Zi(t) at 
each t, a (constrained) Hamiltonian h and a lapse function N(t), it typically has 

*w = *){/w = -|l^ 

SN(t) = -JU(i) (2.2) 

as a gauge symmetry. Here {, } is the canonical Poisson Bracket and e(t) is the gauge 
parameter. 

Because canonical methods will be used, we invoke an abuse of language and re- 
fer to transformations of the form |2.2| as "reparametrizations" throughout this work. 
Reparametrizations of the first type (diffeomorphisms of the t-label set) will be referred 
to as "true reparametrizations." Note that the infinitesimal form of a true reparametriza- 
tion t — > t + St is identical to that of a general reparametrization with parameter e = NSt. 
Because its integrand is a one-form, expression [O] is in fact invariant under the more general 
reparametrizations ( |2.2|) . 

Along a solution on which q = t more than once or not at all, the question "What is the 
value of a(t) when q(t) = r" may be ill-defined. However, expression |2.1| provides a natural 
generalization of the "answer." It adds the values of a(t) at each t such that q(t) = r and 
^ > and subtracts those values of a(t) for which q(t) = r and -4 < 0. Of course, for 



some solutions s the integral in £TT] may not converge, such as when q(t) takes the value r 



an infinite number of times along s, and on such solutions [o]f2 T is undefined. However, this 
will not be a problem for the cases we consider. 

Note that if q — r at several points along a solution, a calculation of [a\ q = T involves 
cancellations between those points with q = r at which dq/dt > and dq/dt < 0. To obtain 
an average without such cancellation while retaining invariance under [2.2| , we might consider 
[b}q=r for 6 = a sign^^) where sign(x) = 6(x) — 6(—x). In particular, [sign^^^] g=r is 
the number of times for which q(t) = r and N(t) > minus the number of times for which 
q(t) = t and N(t) < 0. In a "proper parametrization" for which N(t) > and proper time 
always flows forward, this is just the number of times for which q(t) = r. 

Thus, if we are considering a reparametrization invariant minisuperspace model and q is 
some measure of the size of the universe, the statement that 

lim [signf-J-^l = (2.3) 

pointwise on S (or the on the constraint surface) is equivalent to the statement that such 
models always recollapse; i.e., that any given solution s does not reach arbitrarily large 
values of q. For such a model, we in fact have 

lim [a] 0=T = (2.4) 

T— >00 H 

for any a = a(zi). In [V], we will deduce an analogous property of corresponding quantum 
operators for certain minisuperspace models. We interpret this as a demonstration that our 
quantization of these models describes recollapsing cosmologies. 



III. GENERAL DIRAC METHODS 



We now use the heuristic ideas embodied in |2.1| to quantize time reparametrization 
invariant systems. This section gives the general setting and formalism to be used while 
important examples are described in the next two sections. When the quantum version 
of the integral |2.1| converges, we will define the [^4]q=t through integration and show in 
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! I IB] that they are observables in the canonical sense; that is, that they commute with the 
Hamiltonian. We introduce a space of physical states as in || and find a physical inner 
product in [111 C with respect to which our operators are symmetric. The results of this 
section assume convergence of the integrals, which will be separately discussed for each of 
the examples of [TV] and |V[ However, results for systems outside the scope of [TV] and [V] 
should be considered as formal. 



A. The Setting 

Let us suppose that we wish to quantize a time reparametrization invariant system 
described by a phase space which is a cotangent bundle T*Q over a configuration space 
Q. We will also assume that the dynamics is expressed through a Hamiltonian h and lapse 
function N. Because time reparametrization is to be a gauge symmetry of this system, the 
Hamiltonian is classically constrained to vanish: h = 0. 

If we momentarily ignore the reparametrization invariance and the constraint, this system 
has a straightforward quantization. States of the system live in the Hilbert space L 2 (Q) and 
the Hamiltonian h becomes a Hermitian operator H. In general, we will use capital letters for 
quantum operators and lower case letters for their eigenvalues and classical counterparts. H 
generates a unitary (Newtonian) time evolution and may be used to define a set of Heisenberg 
picture operators corresponding to the functions and vector fields on Q in the usual way. 
In order for H to be self-adjoint, the L 2 (Q) functions may be required to satisfy additional 
boundary conditions. 

Thus, our only tasks are to impose the constraint and to implement the time 
reparametrization invariance. For the second of these, we introduce a lapse operator N(t) 
and an evolution in a parameter time t instead of the Newtonian time evolution described 
above. We proceed with the following algebraic approach: 

For each smooth function q and each vector field p on Q, introduce one parameter families 
Q(t) and P(t) of elements in a noncommutative algebra. Recall that such q and p form an 
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overcomplete set of functions on T* Q and that the classical Poisson Brackets are summarized 
by 

{91,9a} = 
{9i, Pi} = £ Pl 9i = 9'(9i,Pi) 
{PUP2} = {pi,P2}c=p'(PuP2) (3.1) 

where qi, q 2 are smooth functions on Q, pi,p 2 are smooth vector fields on Q, C pi denotes the 
Lie derivative along pi, {,}c is the Lie Bracket of vector fields, and the above expressions 
define the function q' and the vector field p' in terms of q±, pi, and p 2 . We therefore impose 
the relations: 

[Qi{t),Q 2 (t)] = 

[Qi(t),P 1 (t)]=Q'(q 1 ,p 1 ) 

[P 1 (t),P 2 (t))=P'(p 1 ,p 2 ) (3.2) 

on our algebra. We also introduce the family N(t), each member of which commutes with 
every element in the algebra. Further, we require 

-aQ(t)=N(t)[H(t),Q(t)] 



i&P(t) = N(t)[H(t),P{t)] (3.3) 



where H(t) is some symmetric factor ordering of the classical expression for H in terms of 
the g's and p's with these objects replaced by the algebraic elements Q(t) and P(t). While 
this choice of algebra is not unique (see f||U) we will not consider other possibilities here. 

The presence of the lapse makes the dynamics reparametrization invariant. That is, the 
above equations are unchanged by the transformations 

e(t) dQ{t) 



5Q(t) = -ie(t) [H(t),Q(t)] = - 
5P(t) = -ie(t) [H(t),P(t)}-. 



N(t) dt 
e(t) dP(t) 



N(t) dt 

5N(t) = -|e(t) (3.4) 



We may thus consider [3.3| as simply a Heisenberg picture version of the Schrodinger equation 
imposed in ||. 

As in unconstrained quantum mechanics, we now search for irreducible Hilbert space 
representations Ti of |3.2| . For fixed t, they can be constructed in the usual way through 
"coordinate representations" carried by square integrable functions on Q, on which the Q(t) 
act by multiplication by q and the P(t) generate diffeomorphisms along the integral curves 
of p. Since we require H to be Hermitian, the dynamics |3l| identifies representations of 
operators at the times t and t' through 



Q(t) = exp(iHf£ dtN{t))Q{t') exp(-zif /* dtN(tf 
Pit) = exp(iH$ dtN(t))P(t') exp(-iH ft dtN(tj) (3.5) 



where, as usual H = H(t) is actually time independent. Since N(t) commutes with every 
operator, we have N(t) = n(t)l for some real- valued n(t) in any irreducible representation 
for which N(t) is Hermitian. From now on, we will assume that n(t) satisfies 

lim / n(t)dt = ±00. (3.6) 



t'— >±oo Jo 

All such smooth choices of n(t) are equivalent and are related to each other by gauge 
transformations. Note that, given a choice of n(t), this representation is unique to the same 
extent as in usual quantum mechanics. In particular, if the phase space is T*R™ with the 
standard symplectic structure, it is unique by the Stone-von Neumann theorem if weak 
continuity is assumed. 

In order to enforce the constraint in our quantum theory, we will say that "physical" 
states \ip) lie in the physical subspace TC p h ys determined (as in 0) by 

H\iP) = (3.7) 

In general, such states will only be delta-function normalizable. That is, such states will not 
lie in Ti aux at all, but in the dual to the nuclear space associated with the spectrum of H. 
We will deal with this feature later. 
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B. The Quantum Observables 



Our task will be to construct observables; that is, operators which are invariant under |3~^ . 



We will again use operators that correspond to the classical objects [a] q=T defined through: 

[A] Q=T = jT dt ^6(Q(t) - r)A(t) (3.8) 

where, for the moment, we will not worry too much about the factor ordering of the integrand 
or the definition of the delta-function. To show that [^4]Q= r commutes with H, we will prove 
the more general result that time reparametrization invariant operators of the form 

tt = J u(t)dt (3.9) 



for one- forms udt commute with H. We consider only those u>(t) for which |3.9| converges 
when acting on the part of 7i corresponding to some spectral interval [a, b] of H that contains 
the eigenvalue zero. Note that |3.9| does not depend on the choice of lapse function n(t). 
Since udt is a one-form, we assume that it satisfies the equation 



i 



£(^) = N{t)[H(t),!%. (3.10) 
dt^N{t)J Wl w 'iV(t) J 



This is the case when a; is a sum of terms of the form: 

mm^r-^m <s.n) 

and the u^>(t) are built from the Q(t)'s and the P(t)'s so that they satisfy 

-i^uf{t) = N{t)[H{t)^{t)\. (3.12) 

To show that Q commutes with the constraint, we consider a basis of if-eigenstates \E, k) 
with eigenvalue E where k is a discrete label that removes any remaining degeneracy. In 
general, these may not represent normalizable states, but only states that are delta-function 
normalizable. We take (Ei, k\E 2 , k%) = 5^\Ei — E^dfafa where this 8^(E 1 — E 2 ) is to be 
interpreted as a delta-function or a Kronecker delta as appropriate. 
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Consider the matrix elements of jM? between any two states of this type. Note that 
(E^hl^-lE^h) ^exp(i(E 1 -E 2 )^dtn(t))(E l ,k 1 \^r\E 2 ,k 2 ) (3.13) 



N(tY" " -*-V^ 7 Jt> -^'J^^Ntf 
so that we have 

(Eh, h\n\E 2 , k 2 ) = 27i5(E 1 - E 2 )(E U A^-^l^, k 2 ), (3.14) 

using |3~rj| . The delta-function in E\ — E 2 guarantees that our operator preserves the spectral 
subspaces of H, and therefore that 



[Q, H] = when |3.9| converges weakly. 



The reader may be concerned by the presence of the delta-function in E\—E 2 for the case 
where the spectrum of H is discrete. This concern appears to be valid and will be discussed 
in section |VI| Here, however, we will only be interested in applying these techniques to the 
case where the spectrum of H is purely continuous (at least near the zero eigenvalue). 

We now return to observables of the form |3.8| . Our first task will be to address the factor 
ordering of the integrand as j^Q(t) will typically involve P(t) and therefore fail to commute 
with Q(t) and 5(Q(t) — r). From |3.22| and |3.24j , we see that any symmetric ordering of the 
factors -§tQ(t), Q(t), and A(t) will produce a Hermitian gauge invariant that we may take 
to be a "quantization" of [a] 5=T . However, experience shows that the ordering 

[A]Q=r = \ dt{A(t), p(Q(t) - r)} + (3.15) 

is particularly convenient and leads to tractable calculations and analyses. As opposed to 
other factorizations that might be chosen |JL2|| , it has the reassuring property (see [IV B| ) that 
[sign(P )]x°=r is the identity operator for the free relativistic particleQ. Here, {X, Y} + = 
XY + YX is the anticommutator and 9{Q(t) — r) is the projection onto the positive part 
of the spectrum of Q(t) — r. We thus take |3.15| as the definition of [v4]q =t in what follows. 



1 Note that this is the correct result as it is this operator that counts "the number of times that 
c° = t." Thus, [sign(po)L°=r = 1 in the classical case and [l] ;E o =T = sign(po)- 

11 



Before deriving further results, let us reflect for a moment on the two definitions of 
"observable" mentioned above. Within our scheme, the canonical idea that an "observable" 
is an object that commutes with the constraint (the Hamiltonian, H) does not capture the 
entire notion of "gauge invariant" (see related comments in |8||J). For example, the object 
N(t) commutes with the constraint but does not correspond to a gauge- invariant quantity. 
Also, note that the commutator of an operator O with H is not in general gauge invariant 
so that this definition is itself "gauge dependent." For example, an argument similar to that 
of [3.13| and |3.14j shows that any object of the form /f^, dtA(t), where A(t) is an algebraic 



combination of Qi(t) and Pj(t), commutes with the Hamiltonian in a representation for which 
n(t) = 1. However, for more general choices of lapse, it does not commute with H . (This 
problem can be avoided by requiring that [O, H] = in all irreducible representations and 
not just be zero in the representation used.) However, we have seen that reparametrization 
invariant operators do commute with the Hamiltonian. Thus, we take reparametrization 
invariance as the more fundamental notion of observable. 

Finally, we give a few more definitions and useful results. In order to study the case 
where A(t) does not commute with both Q(t) and ^Q{t), we introduce the operators: 

JtA(t)(-6(Q(t) - t)) (3.16) 



and 



Jt{-B{Q{t)-r))A{t) (3.17) 



which satisfy 



= [A] 



R 

Q=T1 



lA] Q =r = l([A] L Q=T + [A}K =T ), (3.18) 

and, if [A(t), p{Q(t) - r)) = 0, we have [A) Q=T = [A} L Q=T = [A]* =T . 
The convenient expressions 

[A] L Q=T \Ek- 1') = 2m J2 \E, k* ; t')(E, k*; t'\A(t') [H, 9(Q(t') - t)]\E, k; t') (3.19) 
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and 



[A]% =T \Ek;t') = 27riJ2\E,k*;t , )(E,k*;t , \[H,9(Q(t') - r)}A(t')\E,k;t') (3.20) 



for the actions of these operators follow from p. 14 



since 



J_ JL 

N dt 



e(Q(t)-T)=i[H,6(Q(t)-r)}. 



Here, we have again assumed that E lies in the continuous spectrum of H. In [IV B| and [V] 
these relations will be used to study L4]q =t . 



C. The Hilbert Space 



Equation |3.14| is also sufficient for us to determine the physical inner product. In the 



spirit of [LJ], we wish to choose this inner product such that all Vt of the form |3J| are 
Hermitian (or at least symmetric) when u(t) is hermitian. Thus, we require 



(\E = 0; h),n\E = 0, h)) P hys = (\E = 0; k 2 ), £l\E = 0, h)) phys 



(3.21) 



where the overline denotes complex conjugation and {.) p hys is the physical inner product. 
But, assuming that the spectrum of H is continuous, we may rewrite |3.14| as 



n\E u h) = ]T \E U k 2 ) (E u feal^Sl^i, h) 



and we see that 3.21 is satisfied if 



\E = 0; k 2 ), \E = 0, h)) phys = 5 kiM 



(3.22) 



(3.23) 



since 



(^ = 0,^1^1^ = 0,^) 



(3.24) 



We conclude that 



All (convergent) observables of the form ^]9| are symmetric on Ti. p h ys with the inner 
product |3~23f| . 



2 This inner product has been independently introduced several times in related contexts. See 
[10,11 for the cases known to the author. 

13 



Uniqueness of this inner product can be guaranteed by enforcing |3.21| for sufficiently 
many Q, such as those (Q± n ) defined by any set of one- forms for which 



(E = 0, k 2 \ N . \E = 0, kx) = (Sk 1 ,j 1 Sk 2 ,j2 ± <W 2 4 2 , 



32 J 



;(l=Fl)/2 



(3.25) 



However, finding a set of physically interesting operators whose Hermiticity makes |3.23 
unique is a more difficult goal which we will not pursue here. 



IV. SIMPLE EXAMPLES 



We now present two simple examples to illustrate how the formalism of |TTT| may be 
applied and to show that it produces reasonable results. These test cases show that the 
methods of |IT| give the usual results for "already deparametrized systems" ( [IV A|) and lead 



to a familiar Hilbert space for the relativistic free particle ( [I V B|) . These results validate the 
methods so that we may then apply them to cosmological models in |V|. 

In the examples we will see that our approach in fact provides an "overcomplete set" of 
quantum observables. By this we mean a set of observables whose classical counterparts are 
overcomplete (see. for example [O) on the (constrained) phase space. 



A. Already Deparametrized Systems 

A standard testing ground for ideas regarding quantization of time reparametrization 
invariant systems is the class of systems that have been deparametrized (i.e., for which a 
time function has been identified Jl|) or, what is nearly equivalent, systems that (in the 
language of canonical quantization) have a Hamiltonian constraint of the form 

Po + h = (4.1) 

where po is some momentum and hi is independent of po. For such systems, the coordinate 
q° conjugate to p effectively acts as a clock and the system may be regarded as "already de- 
parametrized" by the intrinsic time q°. A typical example is a parametrized nonrelativistic 
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particle |[14] i n which case [0] describes such a system with Hamiltonian hi and Newtonian 
time function q°. Note that hi may involve q°, in which case the corresponding Newto- 
nian Hamiltonian is time-dependent. We will see that the approach of [TTI| gives the usual 
quantization of nonrelativistic particles when applied to systems of this form. 

We wish to evaluate matrix elements of [A]qo =t and begin by noting that Q°(t) = 
Q°(t') + T(t) where T(t) = $ N(t) for some fixed t' . Thus, we find 



(<&M-A[A<P=r\<&MJ) = (qi,h;t'\A(t T _ q o)\qlk 2 ;t') 



(4.2) 



where t T _ q o is the time t when T(t) = r — q®. Equation |4.2j can be used directly to show 
that [A]qo =t preserves physical states, but this already follows from the general arguments 
of 

From |4~2], we see that [y4]Q0 =r |^) has components 



T{t)M;t\[A] Q o =T \ 



T(t)M-AA{tm 



(4.3) 



on the Q°(t) = r — T(t) subspace. Now, if \ip) is an eigenstate of H with eigenvalue E, so 
is \ip') = [A]qo =t \iP) and it satisfies 



i^o fei; tW) + h; t\H x {t)W) = E(q°, h; #') 



(4.4) 



It is therefore completely determined by the components given in |4.3| , much as is the case 
in, for example, Since the action of [A]qo =t on \ip) is determined through fO| and p~4] , 
it follows that, for any Z x (t), Z 2 (t) built from the Q(t)'s and P(t)% if [Z x {t), Z 2 (t)] = Z 3 (t) 
then we have 



We can also derive 



[Zi]qo =t , [Z 2 }qo =t = [Z 3 ]qo 



(4.5) 



d r ~ 

1 ^.[ Z Aq°=t = [Hi(t), [Zi}qo =t 



(4.6) 



by considering the matrix elements: 
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(<£ = T-T(t),ki;t\-i-j^[A] ( p =r \i/>) 

= -^(^i =r-T(t),k 1 ;t\[A] Q0=r \iJ>)) - (g? = r-T(t),k 1 ;t\P (t)[A] QO=T \^) (4.7) 

and using [H^t), Q°{t) = 0, [P (t),A(t)] = 0, P Q (t) = Hft-H^t), and [H(t), [A} Q0=T ] = 0. 



In |4.6| , H\{t) is the operator built from the [Q]qo =t 's and the [P]qo =t 's in the same way that 
Hi is built from the Q(t)'s and the P(i)'s. Thus, we have reproduced the usual results at 
the level of algebras. Since the states \E, hi) which satisfy H\E, hi) = E and (q°, k±\E, k 2 ) = 
e tq ° E 5k 1 ,k 2 form a basis for 7i, this also follows at the level of the physical Hilbert space from 
our use of |3.23| . 



B. The Relativistic Free Particle 

Though hardly a sufficient test of a quantization scheme, every procedure for quantizing 
reparametrization invariant systems should be applicable to the relativistic free particle. We 
therefore take this as our next special case and discuss the convergence of the integrals |3.16 



and |3.17| and the operators they define. We will also see that the inner product |3.23| is 



unique, provided that a simple set of observables are required to be symmetric. 

We consider here a relativistic free particle in a flat Minkowski spacetime of an arbi- 
trary number of dimensions. Because the equations of motion can be solved easily, we can 
explicitly write the desired algebra^ as: 



i 



[X' 1 (t).X" (/-')] - f dtN(t), [N(t),X»({)] = (4.8) 



where rj^ u is the Minkowski metric. This follows from |3.3| in the form 



^P»(t) = 0, ^ = ^ (4.9) 
dt M ' ' dt M m " K J 



3 Note that a "deparameterized algebra" based on gauge fixing x° = t may also be introduced for 



this system and can be used in place of 14.81 to construct the operators 3.8. A brief comparison of 



these two possibilities, as well as some of the calculations below, appeared in [12| 
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so that P v (t) = P u . The constraint is given by H = P 2 + m 2 = 0. As in |TJ, iV(t) is some 
function times the identity operator. 

For each t, it is useful to introduce two bases {|:r;i)} and that satisfy 

X»(t)\x;t) =x"\x;t), P u \p;t) = p u \p;t) 

(x; t\x'; t) = 5(x - x'), (p; t\p'; t') = S(p - p') exp(z J dtN(t)p^) 
( X ;t\p;t) = ^=e i v*. (4.10) 

We then proceed to consider the operators [A]xo =T . Inserting complete sets of states and 
using JdEE a \E,a)(E,a\ = J dp\p;f)(p,f\, S(E — E*) = 5(p 2 — (p*) 2 ), and 

(p',t\8(X°(t)-T)\p,t) = 5(p-p^) / _ e - ( Wo) (4.11) 



together with 3.19, we find that 



[A} L x0=T \p;t') = -27TZ J dp' dp*5(p 2 - (p*) 2 )\p*-,t')(p 2 -p'l)(p*-,t'\A(t')\p' ,p;t') 

x l^'^Wa-Po) - -T7 1 d- (4.12) 

2 7r(p' - po) 

Here, p is the collection of spatial components of p so that p = (po,p) and the integration 
over x° has been performed in a distributional sense. Note, again, the presence of the energy 
conserving delta-function which guarantees that our operators commute with the constraint. 
The expression |4.12| contains several singular expressions, all of which are regulated by 



the factor (p'l — pjj) so that we find: 

[A} L x0=T \p;t') = J dp*dp' \p*;t')5(p 2 - (p*) a ))<p*;t|A(t)|i/o l p;t>e* r ^-«'Vo +Po]- (4-13) 
If A commutes with P Q , this simplifies to 

[A\ L x0=T \p;t') = [dW E |q5 , F; t') (W; A~W)\m t')^^ [ep ° + Po] (am) 

J e 6 { + ,-} 2 P0 



where we have introduced the reduced matrix elements {p']t'\A(t')\p;t') p that satisfy: 



(p';t'\A(t')\p;t') = 5(p' -p,W-t'\A{t>)\p;t') PQ (4.15) 
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and where po = yp 2 , — p 2 + p* 2 . If A commutes with Q° : the reduced matrix elements |4.15 
do not depend on p . 



Note that the factor of epo + po in [4.14| minimizes the amount by which this operator 



mixes positive and negative frequency states, though some mixing still occurs. Thus, the 
case of A = Xi does not correspond to a Newton- Wigner operator (which would leave the 
positive and negative frequency subspaces invariant). On the other hand, it may be checked 
that [9(P )signP X l ]xo= T is just the Newton- Wigner operator X % nw (t) associated with the 
coordinate X 1 . 

It now follows that Hermiticity of a few interesting operators guarantees that the inner 
product |3.23| is unique. From |4.14| , we see that 



[P sign(P )]io = >;t , > =po|p;0 (4-16) 

which is hermitian on H so that [P sign(P )]^ 0=r = [-^o si g n (-Po)]xo=r = [-^o sign(P )] x o =T Q. 
If we require [Po sign(P )]xo= r to be Hermitian on the physical space (spanned by |e|p|,p; t')), 
it follows that the positive and negative frequency subspaces are orthogonal. If we also ask 
that [Xj]xo= T and [Pj]x"= T {hj £ {1, 2, 3}) be symmetric, then the inner product is fixed up 
to an overall normalization on the positive frequency subspace and a similar normalization 
on the negative frequency subspace. Additionally, if we require [0]xo= T to be Hermitian, 
where G is the time reversal operator (Q\p;t) = \ — p;t)), these normalizations must be 
identical and the inner product is uniquely determined (up to an overal scale factor) to be: 



(\e\p\,P,t'), \e\p\,p;t')) phys = 2^p 2 + m?5(p - p')8 e>e > (4.17) 



Note that |4.17| agrees with |3.23j since dE = 2\p \dpo. It therefore follows from pi A| that 



all operators of the form [v4]#= T are symmetric (for Hermitian A(t) and B(t)) and we need 
not have chosen x° as the "intrinsic clock." Finally, we verify that we have constructed an 
overcomplete set of observables as densely defined operators on Ti p h ys - This follows since 



4 A similar calculation verifies that [sign(Po)]x°=r = 1 as claimed. 
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4.17| shows that [Xi\xo= T , [Pj}x°=r, [Po sign(P )]x°=T an d [6]xo= T yield normalizable states 
when acting on any e H P h ys for which (|e|p|,p; t'), \ip)) p h ys is a sufficiently smooth and 
rapidly decreasing function of p and thus are densely defined. 



V. RECOLLAPSING DYNAMICS: SEPARABLE SEMI-BOUND MODELS 

The LRS Bianchi IX and Kantowski-Sachs minisuperspaces fall into the class of models 
that we refer to as "separable and semi-bound." These are, roughly speaking, models for 
which the potential is separable and which (classically) describe spacetimes whose homoge- 
neous slices expand, reach a maximum volume, and recontract. They are thus "semi-bound" 
in the sense that they do not classically reach arbitrarily large size. 

Below, we apply the methods of section |T| to such models. As before, we verify the 
convergence ( |V B| ) of the integrals that define our observables and construct a complete set 
of densely defined operators on TC p h ys (|V C|) . In addition, we derive quantum analogues of 
|2.3| and [2.4| , showing that our quantization also predicts recollapsing behavior. 

A. The System and Observables 

Separable semi-bound systems are time-reparametrization invariant models that can be 
described in the canonical formalism with a Hamiltonian constraint of the form h = h^ — hi = 
where ho is of the form p\ + Vo(q°) for some canonically conjugate pair (p , q°) which take 
values in (— oo, oo), h\ does not involve either po or q°, and h\ is a Hamiltonian of the type 
appropriate to a nonrelativistic particle. Thus, we consider the case where the configuration 
space Q is of the form Rx Qi and the phase space is T*Q = T*HxT*Qi. In addition, we will 
ask that V(q°) be smooth, that V (q°) ->• oo as q° -»• +00, V (q°) > 0, and V (q°) < We aqo 
for some positive W and a. This last condition is important only for large negative q°. This 
leads to the following: 
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• Definition: A separable semi-bound system is a Hilbert space H aU x = L 2 (1Z) ® Hi 
together with a self-adjoint operator H = H — Hi such that Hi = 1 eg) for some Hi 
and H = H <g> 11 for /T = + ^(<7o) an d a potential ^(g ) which satisfies (on 

every semi-bounded interval (— oo, g )) the bound V(q°) < We~ aq ° for some W, a > 0, 
but for which V(q°) — > oo as g° — > oo. 

These conditions imply that, when viewed as an operator on the Hilbert space L 2 (R) in a 
"g° - coordinate representation," the spectrum of Hq is purely continuous. The condition that 
^o(<Z°) ~~ * 00 w iH allow us to work with convergent expressions that require no regularization 



or careful treatment. From [[L5[], we see that when described by the appropriate variables 
and with a rescaled Hamiltonian, the LRS Bianchi IX and Kantowski-Sachs models may be 
placed in this form with q° = ln(det g) where g is the 3-metric of a homogeneous slice. 



As described in [III A] , we use the auxiliary Hilbert space L 2 (Q) to carry a representation 



of the algebra defined by |3.2| and |3.3| and require that H = H — Hi be a Hermitian factor 
ordering of h on this space (and similarly for Hq and Hi). From this and the above conditions, 
it follows that H has a purely continuous spectrum with a c (H) = (—00, 00). Note that H is 
time independent and that, since they commute, H and Hi are time independent as well. 

As before, we will make use of a basis in this auxiliary space that is tailored to our 
system. This basis will consist of (delta-function normalizable) states {\E , Ei, k)} that are 
eigenvectors of Hq with eigenvalue Eq and that are eigenvectors of Hi with eigenvalue Ex. 
The (discrete) label k is to remove any remaining degeneracy in our description. Thus, 
H\E ,Ex,k) = (E — Ex)\E ,Ex,k). Note that this set of labels leads to a decomposition 
of our auxiliary Hilbert space as a direct product 7i = Ho <E> TLi through \E ,Ei,k) = 
\E ) <S> \Ex,k) since, for systems of this type, the spectral subspaces of H on L 2 (q°) are 
nondegenerate. 

Also as in ill A, the (constrained) Hamiltonian H and the lapse N(t) define a parameter 



time evolution that, for the canonical variables, is given by 

A(t) = exp(iH I dtN(t))A(t') exp(-iH f dtN(t)). (5.1) 
Jt' Jt' 
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In irreducible representations, N(t) is proportional to the identity: N(t) = n(t)H and we 
assume that |3.6| holds. Also, for each t, our representation has a set of states |g;i) labeled 
by points q in the configuration space Q that define a "configuration representation at time 
f and satisfy (q;t\q';t) = S(q,q'). As in the case of the relativistic free particle, it will be 
convenient to introduce a family of "energy bases" {\E , E ly k; t)} (also labeled by t), where 

\E , Ex, k; t) = (exp[i(£ - Ex) £ dtN(t)}) \E , Ex, k; if) (5.2) 

so that inner products of the form {q;t\E , Ex, k;t) are independent of t. The arbitrary 
phase in the definition of \E , Ex, k) may be used to set \E , Ei,k) = \E , Ei, k; t') for any 
t'. 

Since the spectrum of H includes zero, the set of physical states that are delta-function 
normalizable in the auxiliary inner product and satisfy H\ip) = is not empty. They are of 
the form \E = E\,E\,k] t') where E\ ranges over the spectrum of Hi and k ranges over its 
full set of values. From |3.23|, the physical inner product is just 



(\El, El k* ; t'), \E U E u k- t')) phys = 5 (7) (El E 1 )5 k , ik (5.3) 

where the 5^{E*, E ± ) is to be interpreted as a delta-function or Kronecker delta as appro- 
priate to the spectrum of H\. 

A few brief comments are now in order before deriving the rigorous convergence and 
recollapse results. We again wish to consider the action of the operators [A]s =r . Let us 
suppose that B is, in fact, a configuration variable. Note that 



(E' ,E / 1 ,k';t\6(B(t)-T)\E ,E u k;t)= [ dq(E' , E\, k';t\q;t)(q;t\E , E u k;t) (5.4) 

JQ;b(q)>r 

so that we may write 



(£*, E* l7 k*; t | [A] L B=T \E ,E u k;t) 

= -i27r6{E -Ex- E* + E{) f dE' £ (E -Ex- E' + E\) 

E'i,k' 

x {ElElk*;t\A{t)\E' Q) E'x,k';t) [ (E' , E\, k'; t\q; t)(q; t\E , E l} k; t) (5.5) 
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and, similarly for [A]^ =T . 

Two cases are now of interest. The first occurs when the spectrum of Hi is entirely dis- 
crete. Then, the factor states \Ei) are normalizable and the wavefunctions (E'o, E\, k'; t\q; t) 
fall off rapidly in the directions in which Q is not compact and q° remains constant. Fur- 
thermore, because of the assumption that V (q°) — > oo as q° — > oo, these wavefunctions also 
fall off in the positive q° direction. Thus, if b(q) is any configuration variable for which q° 
is bounded from below in the region of Q with b(q) > r, this factor represents a convergent 
expression. At the other extreme, if the spectrum of Hi is purely continuous (as in [IV B| ), 



then this integral will converge only in a distributional sense, as is appropriate. Note that 
convergence of these integrals guarantees that the matrix element in |5l| (and correspond- 
ingly for [A]qo =t vanishes in the large r limit, verifying the recollapsing behavior of \2.4\ in a 
weak sense. 

B. Rigorous Convergence and Recollapse Results 

We now present two results concerning the convergence of the sums in expressions |3.19 



and [3.20| . These results will show rigorously that a quantum analogue of [2.4| holds. As usual, 
the stronger result follows for the simpler case. 

Let us begin by assuming that B = Q° and that A(t) commutes with both Q°(t) and 
Po(t) so that L4]q 0=t = [A]qo =t = L4]qo =t . Thus, we may again introduce the reduced 
matrix elements {E'i, k'\A\Ei, k) that satisfy 

(E' , E\, k'; t'\A(t)\E , E u k; t') = 6{E' - E ) 1 (E' 1 , k'\A\E u k) 1 (5.6) 



and are independent of Eq. Note that 5.6 defines a "reduced operator" A on the factor space 



Tli. The action of L4]qo =t now takes the form: 

[A} Q0=T \E , Ei, k; t') = -i2n (El - E{)\E -E x + E*, E{, k*; t') x \{E{, k*\A\E u k; )i 

E*,k* 

x / dq (E -Ex + Et\q ) o(qo\E } (5.7) 
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where bras and kets with subscripts and 1 refer to states in TCo and Hi respectively. Note 
that while all terms in this sum are finite it remains to show that the sum itself represents a 
normalizable state in the physical space. Note that this sum will be at best delta-function 
normalizable in the auxiliary space since \Eo,Ei,k) is only delta-function normalizable. 

In order to discuss the convergence of the sum in |5.7| , we would like to find a bound 
M(t,E) such that \I(E,E*,r)\ < M(r, E) where 

I(E,E*,r)= ( (E Q = E*\q )(q°\E = E). (5.8) 

Jq°>T 

We use the fact that the states are delta-function normalized so that I(E,E*,t) = 5(E — 
E*) — J qQ<T {E*\q°) (q°\E) . Because the potential Vo(g°) is exponentially small for large 
negative q°, we approximate the function (q°\E) by plane waves through the expression: 

( q °\E) = -±=(e™° + ase-W) + l-A E (q) (5.9) 

for the appropriate a E with \o.e\ = 1 and some A E which we expect to be small. Here, 
we have set p = \J~E~. Note that we have chosen to normalize the plane waves in a manner 
consistent with the normalization (E\E') = 5(E — E') = ^S(p — p'). 
Thus, we may expand / as 

I = S(E-E*)-—5(p-p*)+ 1 1 



2p An(p + p*) y/pp* 

4:7r(p* — p) y/PP* 

dqo^Sf (5.10) 

qo<r y/pp* 



where the bars denote complex conjugation and 



A E (p*)= o e ? "A B (g u ) (5.11) 

Jq°<T 

The first two terms in |5.10| cancel. 

The bound on the potential may now be used to verify that can be regarded as small. 
By the usual arguments, it follows that A^ satisfies 
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(wT 2+E) ^ = Voiq ° ){qolE) (5 ' 12) 



which has the solution 



A E (q°) = ^e ipq ° f° dq' e- 2vpq '° [*° dq"°e ipq "° (q" \E)V (q"°) (5.13) 

J — oo J—oo 

since, for the right choice of die, A^ satisfies the boundary conditions A# — > and ^jA E — > 



as q° — > — oo. Let us rewrite this expression using |5[ 

A E (q°) = e — / dq ,0 e- 2ipq '° / rfg" e 2 ^"V (g" ) 

\/ A.TT J — oo J— oo 



„„0 „„/0 
47T 



+ ^= f dq ,0 e- 2ipq '° I* dq" a E Vo(q" ) 



OO J — oo 

„0 g/0 



+ e ipq ° T dq ,0 e- 2ipq '° f dq"°e ipq "° A { q"°)V (q"°) (5.14) 

J—oo J—oo 



Let f3\{E) and ^(E) be the maximum values of the first and second terms in |5.14| for 



g° G (— oo,r], which exist since V(<7°) is smooth and bounded on this interval by We aq ° . 
Note that Pi(E), (3 2 {E) — > as — > oo since both describe integrals of fixed L p (— oo,r) 
functions against oscillating exponentials. Here, the integral over q'° is highly oscillatory in 
the definition of ^(E) and the integral over q"° is relevant for (3\{E). 

Further, since (q\E) is smooth and bounded as q — > — oo, |A#| takes some maximal value 



A^ on (-oo, r) by §7Lj and we see that Af < /?i(E) + /3 2 (£) + 7 A^, where 7 = ^e aT is 



independent of £7. Let us suppose for the moment that 7 < 1 so that 



m „ ME) + 82(E) 



Thus, 



\A E (q°) - T° ^°e- 2 ^' f'Vv w "° + a E )V (q"°)\ < Wg+MSt . (5 .l 6 ) 



47T J-00 J-00 ' 1 — 7 

Since /?i + \3i — > 0, this estimate is accurate in the E — > 00 limit and A^(g°) — ► 0. A similar 
argument using |5.14| and |5.15| demonstrates that 



V |A f Lo (5.17) 
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in the E* —>■ oo limit and that the L 2 norm of ^== on (— oo, r) is finite and also vanishes 
for large E*. It follows that I(E, E*, r) — > as E* — > and thus that there exists a bound 
M(E, r) > -E*, r)|. To see this, recall that we know that I(E, E, r) is finite so that we 
may ignore the apparent divergences in |5.10| for p* = ±p. Now, for p ^ p* I is a continuous 



function of E* by |5.10| and f5.13| . The first two terms in |5.10| cancel and the next two terms 



are explicitly decreasing for large p*. The pair of terms involving Ag*(±p) are bounded by 
a multiple of |5.17| and the terms A T E *(±p*) vanish as they are the integral of a smooth L 2 
function against a rapidly oscillating exponential. The last term is bounded by the product 
of and |5.17| . Thus, we have derived the existence of a bound (and in fact of a maximal 
value) for \I(E, E*, r)\, assuming 7 < 1. 

We now consider the more general case and show that, by repeating the approximation 
described by |5.14| , we can follow a similar proof without assuming 7 < 1. Note that if 



we iterate this approximation n times, replacing at each step the A# that appears in the 
integral on the right hand side of |5.14| with the entire right hand side, then the term on the 



right hand side that still contains A# can be bounded by an expression of the form: 

r n W n 

A# / d 2n q n V(q 2j ) < -Jl^e^A* (5.18) 

Jr< qi <q 2 <q :i <..<q 2n ^ (n!)^" 

and so approaches zero as n — > 00. Hence, for any Vq of the form specified, there is some 
finite n such that the bound in |5.18| (which we will call 7„Af / ) is less than Af f so that 
7„ < 1. A bound of the form |5.15| again follows with 7 replaced by 7„ and Pi, 02 replaced 



by the appropriate f3i >n and /?2,n which also vanish in the E — > 00 limit. The rest of the 
argument proceeds as before, and we conclude that I(E, E*, r) is a bounded function of E*; 
that is, 

• There is a function M(E,r) such that M(E,r) > \I(E,E*,t)\. 

From |5.7| and our bound M(E, r), we now see that [A]q =t |.Eo, E\, k; t) is a normalizable 
state for any A such that 

(El - £'i) 2 |(-Ei, k*\A\Ei, k)\ 2 < 00, i.e., such that H^Ei, k)\\ Hl < 00. (5.19) 



E*k* 
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where represents the norm of the state in the factor space H\. This allows us 

the conclude the following result: 



Result 1 : Let A be a densely defined operator of the form ligiAon Ti, aux of a separable 
semi-bound system. Further, let A(t) be the one parameter family of operators defined 



through |5.1| with A(0) = A, and let L4]qo =t be the one parameter family of operators 
on Tiphys defined by [3.18| , |3.19| , and |3.20| . Note that our construction fl3.23|) of 7i p h ys 



identifies H p h ys with 7ii at t — 0. Then, under this identification, if \ip) G Hi lies in 
the domain of [A, Hi], it also lies in the domain of [A]qo =t . In particular, if [A, Hi] is 
densely defined, so is L4]qo =t . 

It now follows that quantization of "separable semi-bound models" captures the classical 
notion of "recollapse" expressed through |2.4| for a(t) independent of q°(t) and Po{t)- As in 
0, this is to be expected since we have used the auxiliary space L 2 (<2). Below, we show 
that the norms of all states of the form [A]qo =t |£'q, Ei, k; t') vanish in the large r limit and 
thus that the operator [v4]qo =t converges strongly to zero. 

Note that we may take the bound M(E, r) on the integral / to be given by the maximum 
value of I(E,E*,t) over E*, which exists by our previous argument. Thus, M(E,t) is 
also given by an expression of the form \J 0>T f(E)\ for / G L 1 (r, oo). Since this integral 
converges, it must vanish in the large r limit and, since the norm of [A]qo =t |£'q, E , a; t') in 
Hphys is bounded by the product of M 2 and the norm given in |5.19| , we find 



• Result 2: For A as in Result 1, [A]qo =t —> as r — ► oo in the sense of convergence on 
a dense set of states; that is, L4]qo =t converges strongly to the zero operator. 

Unfortunately, considering only operators A that commute with Q° and Pq, is not suffi- 
cient to conclusively demonstrate recollapsing behavior. Recall that if a classical spacetime 
s expands smoothly and then recollapses, [a] g o =T (s) represents a difference of a term asso- 
ciated with the expansion and a term associated with the collapse. Thus, our discussion to 
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this point leaves open the possibility that this quantization should be interpreted as describ- 
ing pairs of contracting and expanding spacetimes such that the appropriate terms cancel 
for large r. It would be comforting to find that operators insensitive to the expansion or 
contraction of the universe vanish in the large r limit as well. 

As suggested in section [D], we would like to study the operator [sign(P )]q =T . However, 
for technical reasons we will instead work with 



— e(Q° - T>ig„(P oW Q» - T')—\^ (5.20) 



a 



and derive the analogue of [O] for this operator. Note, first, that the factors of 6(Q — r') 
are completely irrelevant classically for r' < r. Despite this, they will greatly simplify the 
details in what follows. Also, H is time independent and the classical ho is constant along 
a solution. Thus, |5.20| is nearly as satisfactory as [sign(P )]Q =T since, if the corresponding 
classical observable approaches zero along a solution for large r, it tells us that the spacetime 
in question does not reach arbitrarily large values of q°. 

Now, for an operator A (such as |5.20| ) that is constructed only from Q° and Pq, the 
matrix element (Eq, E^, k*;t'\A(t')\Eo, E±, k;t') must vanish unless E* = E\ and k* = k. 
Thus, [A]qo =t |_E , 0) Ei, k; t') must be again proportional to \E , Ex, k; t'). From |3.19| and 13.201 , 
we see (if the spectrum of H\ is discrete) that in fact 

[A] QO=T \E , E u k- 1') = 27r|E , E u k- 1') 

x Im(E , E u k- t'\A(E - H Q )6(Q°(t') - t)\Eq, E u k; t') (5.21) 

where Im denotes the imaginary part. If the spectrum of H\ is continuous, the action 
of [^4]g=r on \Eq, E±, k;t) will not, of course, produce a state that is normalizable in the 
physical inner product. However, essentially the same analysis follows by considering wave 
packets and using the fact that A(E — H )6(Q — r) is of the form O <8> II with respect to 
the decomposition TC = Ho ® 'Hi- 

Thus, our task is to show that this expectation value is finite and vanishes as r — ► oo. 
Note that 9(Q°(t') — t)\Eq, Ex, k; t') is actually normalizable in the auxiliary inner product 
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since (q°\E) decays faster than exponentially for large positive q°. We may express the 
expectation value in [5.21| for the case where A = H ^ +1 9(Q — r')sigii(P )9(Q — t') h ^ +1 as 

-^—mr')\sign(P (t'))e(Q (t') - r )%^|V>(r)> (5.22) 

where \i/j(t)) = 9(Q°(t') - t)\E q , E u k;t') and similarly for \^(r')). Note that is the 



matrix element of a bounded operator between two normalizable states and is therefore finite 
and bounded by some multiple of the norm of \i/;(t)). Again, [A]q =t is densely defined on 
our physical Hilbert space. Also, the norm of decreases to zero as r — > oo, so the action 
of [^4]q= t vanishes in this limit on any finite combination of the states \Eq, Ei, a; t'). Thus 
we find: 



Result 3: O t of 5.20 is densely defined and converges strongly to zero as r — > oo. 



The statement that there may be states to which this does not apply (i.e., that are not 
in the appropriate dense set) is of the same type as the statement that there are states in 
nonrelativistic quantum mechanics with infinite expectation value of the energy. For this 
reason, we feel justified in saying that when Q° has the interpretation of a scale factor 
for a cosmology, our quantization predicts that such spacetimes recollapse. Note that this 
rules out the possibility of an "evolution law" of the type [4.6| . There is no exact unitary 
evolution in r-time for semi-bound separable systems, just as there is no corresponding exact 
Hamiltonian evolution at the classical level. 

C. A Complete Set of Operators 

Our final result will be that a complete set of densely defined quantum observables can 
be constructed by the techniques presented above. In order to prove this general result, we 



consider operators defined somewhat differently than in |3.16| and |3.17| . We will also assume 
that, for each E , Ex, the label k takes only a finite set of values. This is the case for both 
LRS Bianchi IX and the Kantowski-Sachs Model. 
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For any operator B(t) = B(Q(t), P(t)) bounded on H, consider the objects 

[b] q0=t = J ^ dtf{H ) w —e{Q»{t) - r')B(t)e(Q° - T ')jj~lWt d{Q0 - T) ) f{Ho) 

[B] Q0=T = J ^ dtf(H )(-9(Q° - T ))—B(Q°(t) - r')B(t)6(Q ~ O^/W (5.23) 
and note that these satisfy the analogues of |3.18| when B(t) is Hermitian on H. Here, / is an 



as yet unspecified real function that will serve to regulate the above expression. Note that 
the factors of -^±z combine to have no effect in the corresponding classical expression and 
that the ^-functions in Q° — t' are classically irrelevant for r > t' . Thus, the classical objects 
[6] g o =r and [b]q° = r differ only by the factor l/f(h ) 2 . It follows that as B ranges over all 
bounded operators built from the Q(t)'s and P(t)'s, [b] g=T ranges over an overcomplete set 
of functions on the space of solutions. 

The usual calculation shows that convergence of the action of our operator on some given 
state rests in the behavior of the matrix elements 

(\E*, El = Eq, k*), [B] L qo =t \Eo, E x = E , k)) phys = 2nf(E*)f(E ) 

x (E*,El = E*,k*\9(Q - r')B9(Q° - r')-^—9(Q° - r)\E ,E 1 = E ,k) (5.24) 

When the spectrum of Hi is discrete, 6(Q° — t)\E , Ex = E\, a) is some normalizable state 
\ip(Eo, a, t)) in 7i with norm z t (Eq) which is a monotonically decreasing function of r and 
does not depend on a. We will use this function to choose an appropriate /. The case where 
Hi has continuous spectrum can be dealt with in much the same way, using wave packets 
and the factorization 7i = Ho ® Hi, but we will not do so explicitly. 

For the discrete case, fix r' < r and let n(Ei) be an enumeration of the energy levels of 
Hi such that if n = n(Ei), the n-th energy level has eigenvalue E\. Similarly, let m(Ei) be 
the dimension, which we have assumed to be finite, of the subspace of Hi with eigenvalue 
Ex. Now, extend the functions n(E\) and m(Ei) which until now were defined only on the 
spectrum of Hi to all of R, say by defining them to be linear for E% not in the spectrum of 
Hi, and let / be any function such that 

n(E )^/m(E )z T r(E ) 
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It then follows that 

L 4n 2 \ \B\ \ 2 z 2 (E{ 

\(\Eq, E{ = Eq, k*), [B] Q0=T \E , Ei = E , k)) phys \ < 



n(E* )n(Eo)^/m(E )m(E* )z 2 T ,(E ) 

(5.26) 



and 



R Att 2 \ \B\ \ 2 z 2 (E* 

= El } k*) } [B} q0=t \Eq,Ei = E Q ,k)) p h ys \ < 



n{E* )n{E )^m{E )m{E* )z 2 Tl (E* ) 

(5.27) 

where ||5|| is the operator norm of B on 7i. Summing over Eq and k*, it follows that the 

— l — R 

norm of the actions of [B]q 0=t and [B]q 0=t on \E , E\ = E , k) is bounded by 

\\W\E»,E X = E M\l h ys < T^=M^ V-I- (5.28) 

^m(E )n(E )z 2 ,(E )^ ( n *) 2 

and 

|[if !£„,£, = E ,k)]\l h „ < ^^l^J^L- (5 ' 29) 
where n* and E* are related by n* = n(E*). Thus, we have 



• Result 4 Let B be a bounded operator on TC aux of a separable semi-bound system and 
let B{t) be the associated one parameter family of operators on H. aux defined through 
|5.1| with -B(O) = B. Then, for some function / which may depend on the system but 
is independent of the choice of B, the operator 

\B] Q0=r = \{\B] L q°= t + \B]qo =t ) (5-30) 

given by |5.23| defines a bounded operator with domain Tiphys- Furthermore, if B is 
self-adjoint, then so is [B]q 0=t . 

In particular, B may be of the form B = (j4sign(P ) + sign(P )^4) for some bounded A that 
commutes with both P and Qo or we may have B = A for such an A. It follows that we 
have constructed a set of bounded operators [B]q 0=t on H p h ys whose classical counterparts 
[6] g o =r form a complete set of functions on the classical space of solutions (or on the classical 
phase space) of the system. We thus say that 
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The [B]q 0=t form a complete set of bounded quantum observables. 



Note that the bound in |5.26| becomes zero as r — > oo since z T (E ) — > and that the 



bound in |5.27| also vanishes since the sum in |5.27| converges absolutely and each term is a 



monotonically decreasing positive function of r. Thus, we also have 

• Result 5: The [B]q 0=t of Result 5 converge strongly to zero as r — > oo. 

VI. DISCUSSION 

We have seen that the construction of gauge invariants through integration over manifolds 
can be combined with Hamiltonian methods for constrained systems to build a complete 
set of gauge invariant operators that are densely defined in a physical Hilbert space for a 
number of time reparametrization invariant models. A physical inner product was identified 
and this approach was shown to predict the collapse of spacetimes for the appropriate 
Bianchi cosmologies. The primary method discussed involved the use of an auxiliary Hilbert 
space of the form L 2 (<2), that is, containing functions that are square integrable over the 
configuration space Q. 

While the approach can be stated for an arbitrary finite dimensional reparametriza- 
tion invariant model, questions of convergence and spectral analysis have only been ad- 
dressed in narrower contexts. The most general rigorous results were obtained for already 
deparametrized systems ([IV AQ , the free relativistic particle ([IV B|) , and the separable semi- 



bound models of [V|. Results for other systems remain formal, but a forthcoming paper [16 
will address the full Bianchi IX model. 

Although only finite dimensional models were discussed in this work, it might be hoped 
that similar ideas could be applied to full quantum gravity. However, a number of issues 
would have to be faced. These include the choice of representation (which may be constrained 
by the parameterized dynamics of |3.3| in analogy with unconstrained field theory [|17j ,|r8|j ) , 



the choice of algebra itself (see |§[|), and the definition of the products of operators that 
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appear in our integrands (ultraviolet divergences). On the other hand, another interesting 
extension would be to apply these ideas at the (C*) algebraic level without introducing a 
Hilbert space at all. 

Several interesting open problems also remain within the current framework of finite di- 
mensional models. These include the search for a general characterization of the convergence 
of the integrals [3.16| and [3.171 and for the self-adjointness (as opposed to just symmetry) of 



the resulting operators on TC P h ys - It would also be useful to show that such operators form a 
complete set in some inherently quantum mechanical sense (such as that the resulting repre- 
sentation on Hphys is irreducible) and to find a simple set of physically interesting operators 
whose Hermiticity guarantees the uniqueness of the physical inner product. 

Before concluding, a few remarks are in order about the limitations of this approach and 



the delta-functions in p. 14j . The point is that while the technique of introducing an auxiliary 
Hilbert space of the form L 2 (Q) is quite general and can (in principle) be applied to any 
system which has a classical canonical description, technical difficulties may prevent this 
method from being useful. The result may be trivial in two ways. First, it may be that 
when the Hamiltonian constraint is chosen to be a Hermitian operator in the auxiliary space 
its spectrum may not contain the value zero. That is, there may be no solutions to the 
constraint equation that are even delta-function normalizable. The other difficulty is that 
the integral by which we would like to define our gauge invariants may diverge on physical 
states. We have seen that the matrix elements are of the form (E, k\ [A]qo =t \E', k') = 
8{E — E')f(E,E'k,k') for eigenstates \E,k) of H with eigenvalue E and some function / 
so that, if zero lies in the point spectrum of H and a solution to the constraint is actually 
normalizable in the auxiliary space, the action of [A]qo =t on \E,q°) contains a divergent 
factor. While it is possible that this divergence is regulated by the factor /, this seems 
unlikely as we shall shortly see that this divergence has a physical interpretation. 

Note that both of these features tend to occur when the spectrum of H is discrete, such 



as in the models |T^[ where the constraint represents a difference of two harmonic oscillator 
Hamiltonians or in minisuperspace models that recollapse in a finite amount of proper time 
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(/ N(t)dt). This is the case for the LRS Bianchi IX and Kantowski-Sachs models if the usual 
lapse of general relativity is used instead of our "rescaled lapse." A similar feature occurs 
in the recollapsing models of POfl . However, discrete spectra are associated with confining 
potentials as in fT9| , |20|| and, if the dynamics is complete, such systems will generally have 
solutions for which q°(t) = r at an infinite number of times in which case [a] 9 o =r may diverge 
as well. For such systems, one might imagine replacing the integral over t with 

1 /•*+ 

lim —, — : — — r- / dt (6.1) 



and hope to obtain a finite result. Equivalently, we may take 3.22 as the definition of the 



operator Q instead of |3^£ . 

If, however, in a gravitational model, we take as one of our coordinates the parameter 
a = ln(det g), then, because gravity classically admits gravitational collapse, we know 
that the "potential" is not confining in the a — > — oo direction. Thus, at least after a 
proper scaling of the Hamiltonian constraint (in minisuperspace models) we expect that this 
constraint may be turned into a Hermitian operator for which zero lies in the continuous 
spectrum of H in the auxiliary Hilbert space so that the methods of [TTI| will be applicable. 
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